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An instructive method of deriving the matter potentials felt by neutrinos propagating through 
matter on Earth is presented. This paper thoroughly guides the reader through the calculations 
involving the effective weak Hamiltonian for lepton and quark scattering. The matter potentials are 
well-known results since the late 70's, but a detailed and pedagogical calculation of these quantities 
is hard to find. We derive potentials due to charged and neutral current scattering on electrons, 
neutrons and protons. Intended readership is for undergraduates/graduates in the fields of relativis- 
tic quantum mechanics and quantum field theory. In addition to the derivation of the potentials 
for neutrinos, we explicitely study the origin of the reversed sign for potentials in the case of 
. . , antineutrino-scattering. 
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I. INTRODUCTION 



In his famous paper from 1978 1], Wolfenstein showed that neutrino oscillations 0] are altered in the presence of 
matter; a phenomenon which came to be known as the MSW-effect after Mikheyev, Smirnov, and Wolfenstein. This 
occurs as a result of neutrinos experiencing potentials due to charged and neutral current scattering on electrons, 
neutrons, and protons as they propagate through matter on Earth. The explicit expressions for these matter potentials 
can be found in much literature, and a brief sketch for how these quantities are derived is given in which deals with 
■ the special case of charged current electron-neutrino scattering on electrons. Following the outlines of this sketch, this 
paper provides a thorough and pedagogical walk-through of calculating all relevant matter potentials that neutrinos 
experience when propagating through matter in Earth. The calculations include many techniques that are relevant 
when setting out to find expectation values in the second quantized formalism, and should hence be of interest 
for students in the process of being equipped to deal with scattering reactions in the framework of quantum field 
theory. The physics underlying the calculation is also explained to make the derivation more comprehensible, with 
Oh! particular emphasis on how the procedure of finding the potentials can be considerably simplified by observing that 
they are independent of the axial coupling constant, in addition to an analysis of why the potential sign is reversed 
| when considering antincutrinos. Especially the latter argument is of particular value for the intended audience, since 
^*| . practically all literature simply states that the sign is reversed instead of explaining why it is reversed. 



II. NEUTRINO EIGENSTATES AND EFFECTIVE HAMILTONIAN 

The vacuum Hamiltonian describing a propagating neutrino flavor state is as previously mentioned altered in the 
presence of matter. This is readily seen by considering the energy contribution to the Hamiltonian from scattering on 
nuclear components such as electrons and neutrons. First of all, the flavor eigenstates produced in weak interactions 
are coupled to the mass eigenstates through [u e is T ] T = U[f\ v 2 v^Y , where U is the unitary mixing matrix usually 
parametrized as in Ref. |2|, 



u = u 23 u 13 u 12 = 



1 

C 2 3 S23 
-S 2 3 C 2 3 



C13 si 3 e 
1 

i<5cp 



-si 3 e 



Cl3 



Cl2 S12 
-SX2 C12 

1 



(1) 



Here, S; 



and Cij 



cos 9ij and 9ij are the mixing angles between the neutrino states. The 5 C p is the 



CP-violating phase with allowed values S CP € [0,27r]. For S CP S {0,7r} we have CP invariance, while the violation 
is at its largest when 5 CP € {tt/2,3it/2}. For our purposes, <5 CP is set to zero. This is motivated by the fact that 
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CP- violation in the standard model requires a non-zero value of #13, and since experiments have placed stringent 
limits on the value of #13, the effective CP-violation is quite small (see e.g. Ref. Q). Note that we have intentionally 
left out diagonal Majorana phases that are of no consequence for oscillation experiments. 

The evolution of the mass eigenstates ^ M [x) in vacuum is described by 

i±*"(x)=HZ*»(x), (2) 

where Ho — diag(£'i, E 2 , E3). Ei, i = 1,2, 3, are the energy eigenvalues for mass eigenstates. To obtain the Hamilto- 
nian in flavor space, we perform Hq = U H^U^ 1 , where U is the mixing matrix from Eq. When the neutrinos 
propagate through matter, we must add the contributions from scattering on particles present. Accordingly, one 
obtains H = Hq + Hi, written out as 

Hi = H% + H p z + H e z + Hw, (3) 

where the interaction Hamiltonians are given as H z — d\&g(V z ,V z ,V z ), i = n,p,e, and H^y — diag(Vj|,, 0, 0). 
The superscript refers to the scattering component while the subscript indicates which gauge boson mediates the 
reaction, i.e. neutral or charged current. For example, V z represents the effective matter potential due to NC 
(neutral-current) scattering on electrons. Note that we have only included matter potentials arising from reac- 
tions with electrons, protons, and neutrons, since the the concentrations of fj, and t particles is virtually zero on Earth. 

In deriving the matter potentials, we are going to need the effective weak CC (charged-current) and NC interaction 
Hamiltonians 

H-w = -^j=Jwy,Jw-> ( 4 ) 



nf = ^Jpz„. (5) 

The Fermi-constant is Gf = V%g 2 /8m^/, while the currents are defined as 

J z = \ E ^ f 1 ^ 1 - 75) - 2Q, sin 2 0w] *l>u (6) 

i 

where i = (l,i/i,u,d), while if is the accompanying particle isospin and Qi is the particle charge in units of e. Now, 
de=dcos9c + ssmdc- The Cabbibo angle 9c has been experimentally determined to cos^c ~ 0.98, so we shall 
exclude the s quark part from now on and set dg = d. This notation has been adapted from Refs. 0,0. 

The derivation of neutrino matter potentials is certainly not new, and has been studied in e.g. Refs. 0, IE S @- 
This paper contributes to the study of matter potentials by presenting an instructive method to derive the matter 
potentials that applies in the same way to leptons as nucleons. 



III. DERIVATION OF Vfo 

First, we seek the quantity 

v w = (ve(pi,si)e(p 2 ,s 2 )\HwWe(pi,si)e(p2,S2)), (7) 

where Hy/ is the charged-current contribution due to scattering on electrons, omitting the superscript cff for clarity. 
Since we are dealing with elastic scattering, it is fair to assume that the neutrinos and electrons conserve their 
momentum as shown in Eq. Q. The low-energy effective Hamiltonian density relevant for CC v e e scattering is found 
in Eq. Q), namely 

H w (x) = % [e(x)^(l - l5 )v e {x)) [V e {x) lfj {\ - 75 )e(z)] • (8) 
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Nave insertion of this density into 

H w = H w {x)dx, (9) 
Jv 

is not correct. The presence of electrons in a medium leads to two important modifications. First of all, the statistical 
energy distribution of the electrons in the medium is accounted for by integration over the Fermi function f(E e ,T) 
which is normalized to J f(E e ,T)dp e = 1. Secondly, since we do not know the polarization of the electrons, an 
averaging over spins 1/2 ^ s is needed. In total, this corresponds to transforming Eq. (JHJ) to 

x [v e {x)jp(l - 7 5 )e(x)]dp 2 . (10) 
Since only electrons with (p, s) = (j>2 1 s 2 ) will contribute to Eq. (JJJ, we obtain 



[e(x)7 /3 (l - 7 5 )v e (x)] [v e (x)j (l - 7 5 )e(x)] \v e {Pi, s 1 )e{p 2 , s 2 )> 

= 2VE 2 ( P2 ) H2(P2) a ^(P2)^ 2 (p2)7 /3 (l ~ 75K(s)] [ve(xhp(l~l5)u S2 (P2)]We{Pl, s l) e (P2, s 2))- i 11 ) 

Making the identification of the number operator Num S2 (p2) = a\ (p 2 )a S2 (P2), insertion of Eq. (|llfl into Eq. Q 
produces 



v w = (ve{Pi,s 1 )e(p 2 ,s 2 )\ 



G, 



AV2V 



J J f(E e ,T)^^^[u S2 (p^(l~l 5 M 



x [u e (x)j (l - 7 5 )u S2 (p 2 )]dx dp 2 \v e (p 1 ,s 1 )e(p 2 , s 2 )>. 
It should be clear that the v e (x) and V e (x) symbols refer to the second quantized fields 



(12) 



s,p 



2Fw c 



2V"w r 



a s (p)u s (p)e- i ^ + 6t(p) Ws (p)e i ^ 
6 s (p)zJ s (p)e- i ^ + at(p)^(p)e i ^ 



(13) 



where ^> is a lepton or a quark. On the other hand, the state vectors \v e {ph s i))i i = 1,2, are simply representations 
for neutrinos with 4-momentum pi and spin Sj. In order to continue without too many complications, assume that the 
material is isotropic and holds an equal number of electrons with spin up as spin down, i.e. a non-magnetic material. 
This leads to 

Vfr = (v.(pi 1 s 1 )e(p 2 ,s a )&X f J /( ^/ ( )iV ; (P2) Efr-foh^ 1 -75K(*)] 

^ ' S2 

X [F e (x)7/3(1 -7 5 )w S2 (p2)]dx dp2|i/ e (pi,5i)e(p2,52)) 



(i/ e (pi,si)e(p 2 ,S2)|— ^ x 



^e(P2) 



i/ e (a;)7 J a(l - 75)^ (z) 



E[ U ^(P2)7 /3 ( 1 -75))u S2 (P2)]dx dp 2 |^ e (pi,Si)e(p 2 ,S2))- 



(14) 



Here, we have used the Fierz identity to re-arrange the u e and e spinors in such a fashion that makes it possible to 
extract the neutrino-spinor part from the summation over s 2 . The remaining sum is evaluated by 



^[u S! (p 2 )/(l- 75 )» Sl (p 2 )] = Tr{(^ 2 + meh P (l -75)} 

S2 

= P2 Q Tr{ 7 V(l-75)} = 4pf, 



(15) 
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with efficient use of standard trace relations. Eq. (|15f) in synthesis with l|14|l leads to 



G F 

v w = ( v eipi, si)e(p 2 , s 2 )|— 7= x 



f(E e ,T)N e (p 2 ) 



V2 " J J E e (p 2 ) 
x V e (x)j l3 (l - 7 5 )^ e (x)p^dx dp 2 |^ e (pi, s 1 )e(p 2 ,s 2 )} 

a F N e 



= (v e (pi,Sl)e(p2,S2)\ -j= X / V e (x)'Yo(l-j5)l / e(x)dx\Ve(pi,Sl)e(p2,S2)), 



(16) 



where we have exploited the isotropy J p 2 f(E e ,T) dp 2 = and the expression for the total electron density 
J f(E e , T)N e (p 2 ) dp2 = N e . Only integration over x remains, such that we obtain 

v w = ("eipi, si)e(p 2 , s 2 )\ x J V e {x)^ Q {l-^ 5 )v e {x)dx\v e {pi,s\)e(p2,S2)) 



(Ve(pi, Sl)e(>2, S 2 )\ X 



V2 WE Ve 



T r {(^+r7vJ 7 (l- 7 5 )} 



Num sl (pi) 



x ajj (pi)a si (pi) dx|i/ e (pi, si)e(p 2 , s 2 )> 



(t/ e (pi,s 1 )e(p 2 ,s 2 )l x 



4£ , ye dx|i/ e (pi,si)e(p 2 , s 2 )). 



v/2 2V£„ 

Assuming normalized state vectors \v e (p\, si)e(p 2 , s 2 )), Eq. (|17|) reduces to 

Vw- = ^ e x — y dx(i/ e (pi,s 1 )e(p 2 ,s 2 )||i/ e (pi,si)e(p2,S2)) = V2G F N e 



(17) 



(18) 



IV. DERIVATION OF V? 



We now set out to find V% due to v a n, a = e,fj,,T scattering. This reaction is mediated by the Z° boson, so we 
must use the effective Hamiltonian density Eq. JSJ ■ Now, the neutron consists of one u and two d. The u part of the 
Hamiltonian is 



while the d part is 



G f 8 

^=[u(x)7 M (l - 75 - - sin 2 6 w )u(x)] \y a {x)^^{l - 75)^(1)] , 



G 4 

— j= {d(x)^{l - 75 - - sin 2 6 w )d{x)] [^(ajfr^l - 75 W (a?)] • 



(19) 



(20) 



These contributions are to be added in the ratio 1:2 to obtain the relevant Hamiltonian for v e n scattering. In total, 
this gives 



G, 



2V2 



V'„7 M [(l-75- 3 sin 2 6 W 



) - 2(1 - 75 - - sin 2 9 w )]tp n [i/ a (x)7 (tl (l - -f 5 )v a (x) 



G f 

2V2 



{tfnl^i 1 ~ 75)^] [v a (xhf,(l - 75>a( a; )] ■ 



We are left with the effective Hamiltonian 



n z {x) = -^=Jf(E n ,T) x ^^[^(xWil—ftW^x)] [V a (x) 7fl (l-j 5 )v a (x)]dp n . 



(21) 



(22) 



Here, we have introduced the statistical Fermi distribution for neutrons f{E n , T) and summation over the neutron 
spins due to the assumption of unpolarized medium, just as for the electrons. We see that Eq. ill' 21) is of the same 
form as Eq. (|10fl if we use the earlier mentioned Fierz identity, such that the rest of the analysis is equivalent to the 
derivation of Vfy- Following the same procedure as above, we find 
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GfN f 

Vz = -(v a (pi,si)n(p 2 ,s 2 )\ x J V a (x)j (l - j 5 )v a (x)dx\i' a (p 1 , S!)n(p 2 , s 2 )) 

= -{^( Pl , Sl )n(p 2 , S2 )\^- x | Tr{(^ Q +m,J 7 °(l - 7 5 )} 

Num„ 1 (p 1 ) 



x a 



L (PiHi (Pi) dxKOi, si)n(p 2 , sa)> 



(^ a (pi,5i)»(p2,s2)| G f^" x 1 / 4^ l/a dx|i/ (pi,si)n(p 2 ,s 2 )>- (23) 



2^2 2VE U{ 

With normalized state vectors \v a (pi 7 si)n(p 2 , s 2 )) 7 Eq. (|2*5|) leads to 

Vz = - X ^ J &-x{v a (jpi, si)n(p 2 , s 2 )\\v a {pu si)n{p 2 , s 2 )) = - Gp Jl n . (24) 

V. DERIVATION OF VJ 

Now, the expression for V| is found by including the relevant terms from Eq. which for NC v e e scattering reads 

-^f=[v e (x)-f(l ~ 7feK0»0] [e(x)^(l ~ 75 - 4 sin 2 w )e(x)] . (25) 

Taking into account the Fermi distribution of electrons and the averaging of spins, the expression for V§(x) takes the 
form 

Hz ^=-^J /(^e,T)x i^[F e (x) 7 ni~75)^(^][e(x) 7M (l-75~4sin 2 ^)e^)]dp e . (26) 
As before, the rest of the analysis is equivalent to the derivation of V^. In analogy to Eq. I|17|) . it is found that 

Vz = -{i>e(jPi,si)e{p 2 ,s 2 )\ Jr x / (1 - 4sin 2 Bw)v : a (x)j Q (l - 7 5 )^ a (x)dx|i/ e (pi, si)e(p 2 , s 2 )} 



2V2 
G F N e 

2V2 " 2VE Ve 



(v a (pi,s 1 )e(p 2 ,s 2 )\£^ x ^i^- / ( I -4sirT0„ )Tr [(/»,, -:■ ///,, |-."(] 

Num al (pi) 



x (pi)a Sl (pi) dx\v e (pi, si)e(p 2 , s 2 )) 
GpN 1 f 

= -{u e (j) 1 ,s 1 )e(p 2 ,s 2 )\^yJ- x J (l-4sin 2 6»w)4^ e dx|i/ e (pi,si)e(p2,S2)). (27) 

With normalized vectors \v e (px, Sx)n(p 2 , s 2 )), Eq. l(T7)l gives 

11 1 "''11 ) J UA...I'. . .si .( i/'j. l/'x- -N )' W'2- • v ->'/ 

(28) 



^ = — 2~72 =_ X V 7 *" 1 ~ 4Sin ^ 7 ( ^{ l/ e{pi,si)e{p 2 ,s 2 )\\v e (pi 1 si)e(p 2l s 2 )) 
G F (1 -4 sin 2 <V)iV e 



V2 



VI. DERIVATION OF V£ 



In order to derive V% we attack the problem in the same way as for ■ A proton consists of two u and one d quarks, 
so first we find consider the individual contributions from each quark. From Eqs. Ijl9(l and H20JI the total contribution 
to the proton is seen to be 

Gf r 8 4 1 r 

VV7 M [ 2 (1 -75 - g sin 2 9 W ) - (1 -75 - - sin 2 6 w )]ip P V a (x)^(l ~ 7s)^a(^) 

[? p 7 M (l - 75 - 4sin 2 9 w )^ p ] [V a (x) lp (l - 75)^(x)] , (29) 



2^2 
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when added in the ratio 2:1, and including the extra factor of 2 due to the double quark contribution from Eq. |jSJ. 
This gives the effective Hamiltonian for v a p, a = e, /j,,t scattering, namely 

H z( x ) = ^J f(E P ,T)x i^[^(^(l- 75 -4sin 2 ^)^(^)][T7 Q (x) 7A1 (l-75)^W]dp„. (30) 
The further analysis is then just as for Vyy, and leads to 

GjrIV f 

Vf = -{v a {pi,si)p(j>2,s 2 )\ f ^ x / (1 - 4sin 2 6' l y)T7 Q (x)7o(l - 7 5 )^(x)dx|^ Q (p 1 ,s 1 )p(p 2 , s 2 )) 



2^2 

(MPi, si)p{p2,s 2 )\?0?- x J (1 - 4sin 2 9 w )Tr{(^ a + m !/ J 7 °(l - 7 5 )} 



Num sl (pi) 



x al 1 (p 1 )a Sl (pi) dx|i/ a (pi, s 1 )p(p 2 , s 2 )) 
= {v a {Pi,s 1 )p{p 2 ,s 2 )\ x 2 ^ y (l-4sin 2 6» w )4£:^dx|j/ Q (p 1 ,s 1 )p(p 2 ,S2)>. (31) 

As before, we require the state vectors \v a (pi, s\)p{p 2 , s 2 )) to be normalized. Eq. I|31|l then gives 

G N 2 f 
= x —(1 - 4sin 2 9 W ) J dx{u a (pi,Si)p(j) 2 ,S2)\\v a (p-L,Si)p(j)2,S2)) 



V p 



_ G' F ^p(l-4sin 2 9 W ) 

VII. JUSTIFICATION OF TREATING NEUTRONS AND PROTONS AS POINT-PARTICLES 



(32) 



At first sight, it might seem startling that the neutron and proton are treated as point particles with respect to the 
neutrino scattering, instead of taking the well-known nucleon form factors into account. But it turns out that these 
are actually irrelevant as far as the matter potentials are concerned. Notice how all NC potentials are directly pro- 
portional to the corresponding vector coupling constant gv of the scattering particle, i.e. N e oc g v — 2 sin 2 9w — 1/2, 
N n oc g v — —1/2, N p oc g v = 1/2 — 2sin 2 #iy The nucleon vector coupling constants are obtained by quite simply 
adding the quark equivalents for the u and d in the correct ratio; 2:1 for p and 1:2 for n. 

An attempt to apply the same logic to the axial coupling constant <?a would not be successful. The axial current 
anomaly (see e.g. Ref. :9|) modifies the value of this quantity, thus yielding for instance g v ~ 1.27/2 instead of 
1/2 (this particular number has been determined by experiments, see e.g. Ref. One compensates for this 

fact by operating with form-factors in scattering reactions on nucleons. For neutrino forward scattering (essentially 
the MSW-effect), one needs to include form factors that are non- vanishing in the limit of zero momentum transfer. 
Because of the mentioned axial current anomaly, the form factor for the axial current will give results different from 
a combination of quark currents in the respective ratios for neutrons and protons. But the potentials were generally 
derived to be independent of gA, regardless of scattering component. As a consequence, we need not worry about the 
axial form-factors. Left with only the vector part, the nucleon can simply be treated as the superposition of three 
quarks in this scenario. Of course, such an uncomplicated picture is certainly the exception rather than the rule with 
respect to nucleon scattering. 

VIII. HOW TO TREAT A VARYING DENSITY 

Eqs. I|24|) . I|28|) . ll-il'l) are equally valid even if the particle density Ni fluctuates in space. One would compensate for 
such behavior by evaluating the matter potential at each separate point in an artifical lattice-space, or alternatively 
approximating the particle density with a contineous function of the space-coordinate x. The two cases constant vs. 
non-constant particle density differ crucially in terms of how they affect the probability for a neutrino oscillation to 
have occured (see e.g. Ref. for an introduction to the topic of neutrino oscillations). Dependent on which of the 
two scenarios that is applicable, the oscillation probability takes on a certain characteristic spatial-dependence. This 
is also dependent on the nature of the particle density fluctuations, i.e. adiabiatic, unit-step profile, et.c. However, 
oscillations are not the topic of this paper, and we shall settle with the derivation of Eqs. Q24[l. (|28|l . (|32|l since these 
are heuristically valid even in the case of a varying particle density. 
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IX. POTENTIAL SIGN FOR ANTINEUTRINOS 



The matter potentials derived so far has been for neutrinos. One cannot generalize these results to be valid for 
antineutrinos, and here is why. Consider the induced matter potential for v e due to NC scattering on electrons. This 
potential is given as (see Sec. 01 



Vz = (f e (pi,si)e(p 2 , s 2 )\H e z \v e {pi, s 1 )e(p 2 ,s 2 )), 
while the induced matter potential for V e for the same type of reaction reads 

V% = (J 7 e(Pi,si)e(p2,s 2 )|F||Fe0i,si)e0 2 ,s 2 ))- 
In the general case, the operator part of Eq. 1|34[1 takes the form 

(o|Mp)[EMk)&t(k')]^( P ')|o), 

k,k' 



(33) 
(34) 

(35) 



where the sum over 6„(k)&|,(k') comes from the second quantized fields v e (x)v e {x). Using the commutation relation 
[b v (k),bi(W)] = 5(k' - k), Eq. g3 becomes 



-(o\b„( P )[Y,bi(k>)b„(k)]bi(p')\o). 



(36) 



k,k' 



The only non- vanishing contribution from this term is obtained by taking k = p',k' = p. In our case of elastic 
scattering, we impose the limit p — * p'. A similar argument for Eq. I|33[) produces the operator sequence 



Mp)[E4(kOa„(k)]at(p')|0>. 



(37) 



k,k' 



Here, a„(p) and aj,(p) (b v (p) and &£(p)) are interpreted as annihilation and creation operators for neutrinos 
(antineutrinos). The matter potentials differ in sign due to the relative minus sign between Eqs. (|3"F))l and (|3"7j) . Thus, 



V z = ~V§. 



X. SUMMARY 



The relevant matter potentials for neutrinos propagating through Earth have been derived using a method that applies 
in the same way to both leptons and nucleons. We have explicitely shown why the matter potential sign is reversed in 
the case of antineutrinos. Specifically, it is worth to note that the effective potential in an electrically neutral medium 
with respect to neutrino oscillations is given by in Eq. (J3J. From Eqs. 128|) and 132fl it is easily seen that 



vz + v* 



= 0. (38) 

N e =N„ 



Thus, there is no effective matter potential felt by neutrinos due to scattering on protons and electrons mediated by 
Z° in an electrically neutral medium. Also, since H z in Eq. (|3J) is diagonal, it gives an overall phase shift to all 
neutrino flavors, which is of no significance in the oscillation scenario. Our results are summarized in Tab. [I] 



TABLE I: Neutrino matter potentials induced by Earth. 



Type of reaction 


Matter potential 


Vz 


TG F N n /V2 


vi 


±Gf (1 -4 sin 2 6w)N p /V2 


UJ 


TG f (1 -4 sin 2 6 w )N e /V2 




±V2G F N e 
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Generalization to matter potentials such as V% are obtained simply by substituting N e — > N^. The upper sign refer 
to neutrinos, while the lower gives the matter potential for antincutrinos. 
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